The Chiral Magnetic Effect 
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Topological charge changing transitions can induce chirality in the quark-gluon plasma by the 
axial anomaly. We study the equilibrium response of the quark-gluon plasma in such a situation to 
an external magnetic field. To mimic the effect of the topological charge changing transitions we will 
introduce a chiral chemical potential. We will show that an electromagnetic current is generated 
along the magnetic field. This is the Chiral Magnetic Effect. We compute the magnitude of this 
current as a function of magnetic field, chirality, temperature, and baryon chemical potential. 



I. INTRODUCTION 

The quark-gluon plasma is a phase of extremely hot 
matter consisting of quarks and gluons. Just after the 
Big-Bang, the universe itself was in the quark-gluon 
plasma phase. The quark-gluon plasma can be created 
and studied using collisions of heavy ions. An active ex- 
perimental program to investigate the properties of this 
hot phase of matter is underway using the Relativistic 
Heavy Ion Collider (RHIC) at BNL. In the near future 
the quark-gluon plasma will also be studied using the 
Large Hadron CoUider (LHC) at CERN, the FaciHty for 
Antiproton and Ion Research (FAIR) at GSI, and the 
NICA faciUty at JINR, Dubna. 

The behavior of the quark-gluon plasma is described by 
Quantum Chromodynamics (QCD). One of the intrigu- 
ing predictions of QCD is that in the quark-gluon plasma 
phase certain special gluon configurations to which one 
can assign a winding number play a role ^, '2'|. This 
winding number is a topological invariant, which means 
that smooth deformations of these configurations do not 
change the winding number. Experimental evidence for 
the existence of configurations with nonzero winding 
number is only indirect from the meson spectrum [3|,[3,0]- 

The configurations with nonzero winding number are 
in fact transitions which invoke passing a potential bar- 
rier with a height of order the QCD scale Aqcd over the 
strong coupling constant as- Because of the height of the 
barrier, the transitions are highly suppressed at low tem- 
peratures since they require tunneling 0. The configu- 
rations responsible for this tunneling process are called 
instantons [3,[3,[^,0|- At high temperatures in the quark- 
gluon plasma phase, it is possible to jump over the poten- 
tial barrier. The transitions are therefore not su ppr essed 
anymore and called sphalerons [1, H, [l^, [ll|, H^ . [isj . 
These configurations were studied in the electroweak the- 
ory as a mechanism for baryogenesis [13, [HI El, , and 
are also relevant for QCD [l5l [TgI. [l7|. 



At these high temperatures the configurations with 
nonzero winding number can be produced with relatively 
high probability [HI, [ll|. Therefore the quark-gluon 
plasma is the best place to find direct experimental evi- 
dence for the existence of gauge field configurations with 
nonzero winding number. 

These configurations do something very distinct to 
quarks; they can, depending on the sign of their winding 
number, transform left- into right-handed quarks or vice- 
versa via the axial anomaly [19,] (see also [H, HH). For 
massless quarks, the axial anomaly equates df^j^ to the 
topological term. The spatial integration of cJ^j's yields 
an exact relation for the rate of the chirality change in- 
duced by topological configurations, which reads 
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where Nh^l denotes the net number of quarks (minus an- 
tiquarks) with right/left-handed chirality, Nf the num- 
ber of massless flavors, and F^j, — \tfiu\aF^'^°' , with 
g0i23 ^ 2^ i-jjg massless flavors equally couple to the 
gauge field, hence the proportionality factor Nf arises 
in Eq. ([T]). Let us stress that, in the common conven- 
tion, chiral quarks have opposite helicity to antiquarks; 
a particle with right-handed chirality has right-handed 
helicity, while an anti-particle with right-handed chiral- 
ity has left-handed helicity. For instance the helicity of 
the antineutrino is right-handed. Here right-handed 
helicity means spin and momentum parallel, while left- 
handed helicity means spin and momentum anti-parallel. 
Therefore the difference — Nl can also be read as the 
total number of quarks plus antiquarks with right-handed 
helicity minus the total number of quarks plus antiquarks 
with left-handed helicity. For physical gluon configura- 
tions (configurations with finite action) the time integral 
over the right-hand side of Eq. ([1]) is equal to minus twice 
the winding number of the gluon field configuration. As a 
result of the axial anomaly the interactions between these 
configurations and the quarks break the parity {V) and 
charge-parity {CV) symmetry. Ordinary (perturbative) 
interactions between quarks and gluons cannot induce a 
difference between the number of right- and left-handed 
quarks. A mass term always will tend to wash out such 
difference 12211. 
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In QCD, the probability to generate either a gluon 
configuration with positive or negative winding number 
is equal. This is because there is no direct V and CV 
violation in QCD (assuming the value of the 6 angle 
is equal to zero). In the quark-gluon plasma, many of 
these configurations can be generated at different points 
in space and time with different winding numbers. In 
pure SU(iV) Yang-Mills theory this process is completely 
random; the dynamics of the chirality change is that of 
a one-dimensional random walk. In QCD with massless 
flavors, however, it will cost energy to induce a differ- 
ence between the number of right- and left-handed quarks 
because of the Fermi-principle. Therefore the dynamics 
is not completely random anymore, and there is a pref- 
erence to decrease the chirality '15]. In any case, the 
variance of the chirality will be nonzero, and increase as 
a function of time according to diffusion. Hence, it is 
expected that every time the quark-gluon plasma is pro- 
duced, it will posses a non-zero chirality 23]. The chi- 
rality averaged over many events of quark-gluon plasma 
production vanishes. Therefore one speaks in this case of 
event-by-event V- and CT'- violation. 

Next to the sphaleron transitions, chirality could also 
be introduced in the quark-gluon plasma in the same 
way due to chromoelectric and chromomagnetic fields in 
the initial state of the quark-gluon plasma produced in 
heavy-ion coUisons, i.e. the so-called glasma [H, [H, |2^. 
Although the net topological char ge c annot develop with 
the Boost invariant configuration [23], the glasma insta- 
bility spontaneously breaks the Boost invariance 12Q], so 
that the event-by-event topological charge fluctuation is 
expected. The situation is then quite reminiscent of the 
sphaleron transitions [2^ . 

It was first argued by one of us 27] that if V- and 
CP-violating processes are taking place in the quark- 
gluon plasma produced in heavy-ion collisions, then pos- 
itive charges should separate from negative charges along 
the direction of angular momentum of the collision. In 
Ref. { 2^ this mechanism was worked out in more detail 
using an effective 9 angle to mimic V- and CP-violating 
processes. In heavy-ion collisions the magnetic field is 
pointing in the direction of angular momentum. It was 
shown in Ref. [2^ that this magnetic field induces charge 
on a 0-domain-wall in such a way that an electric field is 
created parallel to the magnetic field. In this way posi- 
tive charge is separated from negative charge along the 
magnetic field. In Ref. [29] a different mechanism for 
charge-separation was discussed (see also [13]). It was 
shown that a magnetic field in the presence of imbal- 
anced chirality induces a current along the magnetic field. 
Again, as a result, positive charge is separated from neg- 
ative charge along the magnetic field. This is called the 
"Chiral Magnetic Effect" . 

Due to the separation of charge along the direction of 
the magnetic field in heavy-ion collisions, an asymmetry 
between the amount of positive/negative char ge a bove 
and below the reaction plane is expected [13, IH, |2^. 
These asymmetries can be analyzed in experiments using 



a correlation study as proposed by Voloshin 'si'] . Prelim- 
inary data from the STAR collaboration is presented in 
Refs. [13, [111 . Observation of the Chiral Magnetic Effect 
will be direct experimental evidence for the existence of 
topologically non-trivial gluon configurations. It further- 
more is evidence for event- by- event P- and CP-violation. 

The Chiral Magnetic Effect could be used to deter- 
mine whether a deconfined chirally symmetric phase of 
matter is created in heavy-ion collisions [2^. Deconfine- 
ment is a necessary requirement for the Chiral Magnetic 
Effect to work, since it requires that soft quarks can sep- 
arate over distances much greater than the radius of the 
nucleon. Moreover, chiral symmetry restoration is essen- 
tial, because a chiral condensate will tend to erase any 
asymmetry between the number of right- and left-handed 
fermions. 

In this article we will investigate the Chiral Magnetic 
Effect in detail. In order to treat the non-vanishing chi- 
rality, we introduce a chiral chemical potential, denoted 
as ^5 . The chiral chemical potential will be generated by 
the topological charge changing transitions. We will not 
study this dynamical process, but we just assume this 
chemical potential is there. Then we will study the im- 
plications of applying a magnetic field to a system with 
nonzero chiral chemical potential in equilibrium. We will 
see that an electromagnetic current will be induced in 
the direction of the magnetic field. We will compute the 
magnitude of this current as a function of magnetic field, 
chirality, temperature, and baryon chemical potential. 

Besides the Chiral Magnetic Effect, the fact that 
a magnetic field can infiuence QCD processes is well 
known. For example a magnetic field can induce chiral 
symmetry breaking |34j . influence the chiral condensate 
[Hidl], and therefore modify the phase diagram of QCD 
(see Refs. [13, HI] for recent discussions). Also the color- 
superconducting phases predicted to exist at high baryon 
densities are strongly affected by a strong magnetic field 
[39I . I40I . l4ll . [4^ . Finally the anomaly in the presence of a 
magnetic field can give rise to all kinds of interes ting ef- 
fects, like spontaneous creation of axial currents [43114^ 
and formation of Tr^-domain walls [Ssj . 

The analysis we present in this article can be used to 
make predictions for the charge asymmetries in heavy-ion 
collisions like are done in Ref. [29] . We will encounter the 
beautiful physics of the anomaly, current quantization 
and the index theorem, and periodic oscillatory behavior 
due to Landau level quantization. 



II. CHIRAL CHEMICAL POTENTIAL 

As was argued in the introduction, topological charge 
changing transitions can induce an asymmetry between 
the number of right- and left-handed quarks due to the 
axial anomaly. In order to study the effect of this asym- 
metry we introduce a chiral chemical potential This 
chemical potential couples to the difference between the 
number of right- and left-handed fermions. To the La- 
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FIG. 1: Spectrum of massless Dirac fermions with right- and 
left-handed chirahty in the presence of an chiral chemical po- 
tential /is. The subscript ± denotes the eigenvalue of the 
spin in the z-direction. The chiral chemical potential induces 
a nonzero density of right-handed particles and left-handed 
anti-particles. 



grangian density the following term is added 

/i5Vi7°7'^- (2) 

The energy spectrum of the free Dirac equation in the 
presence of a chiral chemical potential is for massless 
modes (with p^: ~ py ~ for simplicity), 



WL± = TPs + M5- 



(3) 
(4) 



Here ± represents the spin in the z-direction and R, L the 
chirality. The momentum in the z-direction is given by 
Ps; let us stress that in our notation ps does not denote 
the third component of a four- vector with metric conven- 
tion. We have displayed the massless energy spectrum in 
Fig. m In the massless limit one can distinguish modes 
with right-handed chirality from modes with left-handed 
chirality. It should be mentioned that p^ is restricted to 
be positive for the R+ and particle modes so that 
the helicity is positive for R+ and negative for L— , re- 
spectively, and p3 is negative for the R— and L+ particle 
modes (see Fig. [1]). If the chiral chemical potential is 
positive some of the right-handed particle modes will be- 
come occupied while some of the left-handed anti-particle 
modes will be filled as well. A net chirality is created in 
this way. 

The chiral chemical potential lifts the degeneracy be- 
tween modes with right- and left-handed chirality. A dif- 
ference between the total number of particles plus anti- 
particles with right-handed and left-handed helicity is 
created. The magnetic field will lift the degeneracy in 
spin depending on the charge of the particle. Hence par- 
ticles with right-handed helicity will tend to move op- 
posite to anti-particles with right-handed helicity. As a 
result an electromagnetic current is generated along the 
magnetic field, which is the Chiral Magnetic Effect [1^ 
(see also Refs. t29|,|30|] for a pictorial representation of the 



Chiral Magnetic Effect). We will compute this induced 
electromagnetic current in the next section. 

The effect of a finite amount of topological charge 
change can also be mimicked by an effective theta an- 

t which could depend on space-time (see for example 
.mm)- One adds to the Lagrangian of QCD the 
following term, 
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(5) 



327r2 

By performing an axial U(l) rotation this term can be 
transformed into the following fermionic contribution 

1 



2Nf 



(6) 



Identifying this with Eq. Q we see that /i5 — do9/2Nf. 
We can also identify fi^ with the time component of an 
axial vector field A^^. The effective theta angle results 
in a difference between the rates of changing left-handed 
into right-handed and changing right-handed into left- 
handed particles. The chiral chemical potential, how- 
ever, is a more static quantity; it is the energy neces- 
sary to put a right-handed quark on its Fermi surface 
or to remove a left-handed quark from its Fermi surface. 
It describes the difference between the number of right- 
and left-handed fermions. An effective theta angle to 
describe spontaneous V and CP-violating processes has 
been discussed often in the literature (for examples see 
Refs. [13, m, H mi)- The chiral chemical potential 
has on the other hand only been used in a few papers 

Let us finally point out that the chiral chemical poten- 
tial has no sign problem, i.e. the fermionic determinant 
with /i5 is real and positive. In the presence of a chi- 
ral chemical potential the fermionic determinant reads in 
Euclidean space-time, 

detA4(M5) =det(|)-|-^57^7^ + m) , (7) 

where p = j'^D^. Here we have chosen a representation 
in which all "je matrices are Hermitian, "Je — l^ilh ~ 
ij^ . Since p and 7^7^ are anti- Hermitian the eigenvalues 
of A^(/i5) are of the form iXn+m, where A„ € R. Because 
75 anticommutes with p + /i57^7^, all eigenvalues come 
in pairs, which means that if jA„ -I- m is an eigenvalue, 
also —iXn -l- TO is an eigenvalue. Since the determinant 
is the product of all eigenvalues we see that the deter- 



minant is the product over all n of A„ 



Hence the 



determinant is real and also positive semi-definite. This 
is very interesting because it allows for a lattice QCD 
simulation of chirally asymmetric systems. The lattice 
QCD can then simulate the Chiral Magnetic Effect by 
introducing a space-dependent phase on the link variable 
which amounts to the external magnetic field. 



III. COMPUTATION OF INDUCED CURRENT 

In this section we will show if a magnetic field is ap- 
plied to a system with an asymmetry between the number 
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of right- and left-handed fermions, an electric current is 
induced along the magnetic field. We will compute this 
current in four different ways, since we think they are 
all very instructive. The first way is through an energy 
balance argument. Then we will arrive at the result by 
solving the Dirac equation. The third way is by explicitly 
computing the thermodynamic potential in the presence 
of a magnetic field. The last derivation we discuss is us- 
ing a derivative expansion of the effective action. We will 
compute the current for a generic fermion with electric 
charge e and neither flavor nor color; at the end of this 
section we will discuss what happens if we recover the 
flavor and color for quarks. 

Let us set up notation. We will take the metric g^^ — 
diag(l, —1, —1, —1) and the chiral representation for the 
gamma matrices; 



7^ = 



cr'" 





7^ = 



-1 
1 



(8) 



where = (1, tr*) and — (1, — cr*) arc the quaternion 
bases. Using this convention it is possible to write the 
fermion field ip into its left- and right-handed components 
"0 = {'Pl,4'r)'^- We define the right- and left-handed 
chemical potentials as = h + fJ,5 and hl = ^ Ms- 
Here /i denotes the quark chemical potential, which for 
three colors {Nc = 3) is equal to one third of the baryon 
chemical potential. If we write ps, as we mentioned in 
the previous section, we mean the z-component of the 
momentum vector p and not the third component of the 
four- vector p^ . 

The total current is equal to the volume integral over 
the current density. 



(9) 



Let us consider fermions with positive charge e in a 
background magnetic field B. The fermions will occupy 
Landau levels, so their motion in the transverse (to the 
field B) plane will be restricted. The fermions however 
are free to move along or opposite the direction of B; 
since the spins of the fermions are preferentially aligned 
along the field, the motion parallel to B corresponds to 
the right-handed fermions, and anti-parallel to S - to 
left-handed fermions. 

The presence of an electric field E parallel to B causes 
the chirality to change (see Ref. [13] for a related discus- 
sion of particle acceleration in cosmic strings). The en- 
ergy of right-handed fermions moving along the electric 
field under the influence of the Lorentz force will grow 
linearly with time, leading to the growing Fermi momen- 
tum. 



p§ = eEt. 



(12) 



Likewise, for left-handed charges the Fermi momentum 
will decrease, with pp = —pp; this corresponds to the 
production of left-handed anti-particles with charge — e. 
Therefore the particles with charge e will move along the 
field, and anti-particles with charge — e, against the field. 
Thus, an electric current is created along E. 

The density of right-handed fermion states is equal 
to the product of the longitudinal phase space density 
dn/dz = Pp/2tt and the density of Landau levels in the 
transverse direction d^n/dxdy — eBjl-n^ 



2tt 



eB 
2^ 



47r2 



E B t. 



(13) 



The same expression yields also the density of left-handed 
anti-fermion states; therefore, the rate of chirality A^s — 
Nii — Nl generation per unit volume per unit time is 
then given by 



The current density is given by the following expectation 
value. 



j^'ix) = e{^{x)^^'i^{x)). 



(10) 



Here the expectation value is over a thermodynamic en- 
semble. One can write the current density in terms of 
right- and left-handed spinors as 

f{x) = e{4>U^)a^Mx)) + e{4>{{x)d^^cl^L{x)). (11) 



A. Axial anomaly and the energy balance 

The easiest way to obtain the right expression for the 
current is using a beautiful argument of energy balance 
by Nielsen and Ninomiya [i^ . Consider a situation with 
an electric field E and a magnetic field B in the presence 
of a chiral chemical potential. In that case the electro- 
magnetic anomaly will tell us that the rate of change of 
chirality is equal to the volume integral over e^E-B /2tt^ . 
There exists an intuitive derivation of this rate [46] which 
we will repeat here. 



We have thus reproduced the general anomaly relation 
for the electromagnetic fields. 

Consider now the energy balance related to the chi- 
rality change. To change a left-handed fermion in a 
right-handed fermion requires removing a particle from 
the left-handed Fermi surface and adding it to the right- 
handed Fermi-surface. This will cost an energy /Xij—^Xi = 
2/^5 or /15 dA^5 . So multiplying this energy by the rate of 
chirality change we know how much energy is needed per 
unit of time. This energy has to come from somewhere, 
assuming no losses; it will be equal to the power delivered 
by a current. This power is equal to the product of the 
current with the electric field. So one finds |49| 



d^xj ■ E = 



dN5 



27r2 



d-^xE-B. 



(15) 



We can take E in the direction of B in this expression. 
Then if we take the limit E ^ wc find 



J = 



27r2 



d^a;^. 



(16) 
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This derivation clearly shows that not only the axial 
anomaly of QCD plays a role in the Chiral Magnetic 
Effect, but also the electromagnetic axial anomaly. The 
QCD anomaly provides the chirality, the electromagnetic 
anomaly the current. In a box with periodic boundary 
conditions, the number of Landau levels is an integer. 
This gives rise to current quantization as we will closely 
see in the next microscopic derivation. 



B. Dirac equation 

We will now compute the induced current by solving 
the Dirac equation in the presence of a magnetic field and 
chiral chemical potential. We take the magnetic field in 
the z-dircction, 



(17) 



The Dirac equation in this background reads 

(i7^P^ - m + Ai7° + Ai57%') ^{x) = 0. (18) 

where D^^ = dp, — ieAp,. In order to incorporate the 
magnetic field given in Eq. P7| the only non-vanishing 
components of are /i = 1,2. Furthermore only 
will depend on x and y. The precise form of A^ is not 
relevant for our calculation. 

We will compute the total current in the z-direction as 
is given in Eq. ^ starting from Eq. (fTTjl . To proceed one 
has to make a momentum decomposition of the fields in 
terms of creation and annihilation operators. As is shown 
explicitly in Ref. [l^l the only non- vanishing contribution 
to 



(19) 



arises from the transverse zero modes, i.e. modes which 
have Px = Py = 0- The reason is that in all the non-zero 
modes there is a spin degeneracy in energy, which results 
in the cancellation of the expectation value of as [3, [s^] . 
The transverse zero-modes are however not degenerate. 
Let us denote the number of transverse zero modes with 
fJa equal to ± as N±. One shows that the difference 
— is equal to the index of a two-dimensional Dirac 
Hamiltonian in the presence of a magnetic field [H^l • This 
index can be expressed in terms of the total flux One 
finds S, [Hi 



^4 



iV_ = 



.27r . 



(20) 



where we have introduced the floor function [a;J which 
is the largest integer smaller than x. The flux is equal 
to the integral of the magnetic field over the transverse 
plane. 



/ 



$ = / d'^xB{x,y). 



(21) 



Let us stress here that the number of zero modes is quan- 
tized, and not the magnetic flux itself. 

It is now possible to construct the total current. It is 
equal to the sum of number densities in the transverse 
zero-mode weighted by the spin degeneracy ±iV-|-. For 
the right-handed modes we find 



nip3 + Hr) 



dp3 
2^ 



n{p3 - 



2n 



(22) 



Here denotes the length of the system in the z- 
direction and n{uj) — [exp(ci;/T) -I- 1]~^ is the Fcrmi- 
Dirac distribution. The two Fermi-Dirac distributions 
in Eq. (|22p correspond to right-handed particle and an- 
tiparticle modes respectively. In front of the antiparticle 
contribution there is a minus sign, since (p^cj) is the num- 
ber density of particles minus antiparticles. The temper- 
ature dependence has dropped out from Eq. (j22p without 
approximation. The reason why p^ runs only positive is, 
as we have explained on Fig. [U R+ has positive ps only 
and R— has negative p^ whose sign we changed in the 
integral. Similarly, for the left-handed modes we find, 



n{P3 + y^h) 



±N. 



n{p3 - (il) 



2tt 



(23) 



By taking the spin sum and subtracting L from R con- 
tributions we find that the total current becomes 



J = e 



27rJ TT 



(24) 



The result is independent of temperature and density. By 
adding the two contributions one finds the total induced 
axial current, J5 = J d^x in the massless limit, 



e<i> 
2^ 



TT 



(25) 



This current was computed for /is = by Metlitski and 
Zhitnitsky We recover the result of Ref. and 

find that the total axial current is independent of /X5. 

This derivation can be performed in the more general 
case with massive fermions. The computation is more 
involved, but the final answer will turn out to be inde- 
pendent of mass. In the next derivation we will include 
a mass term and show that the answer is independent 
of mass. There, it will be transparent why the result is 
insensitive to temperature and density as well. The last 
derivation using the derivative expansion will provide un- 
derstanding why the current is independent of mass from 
a different point of view. 



C. Thermodynamic potential 

We will now derive the current in a homogeneous mag- 
netic background using the thermodynamic potential. In 
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the presence of a chiral chemical potential we find that 
the thermodynamic potential is given by 



\eB\ ^ ^ 

s=± r!=0 



dp3 
27r 



+ r^log(l + e-'3('^''-±^); 



(26) 



where n is a sum over Landau levels, s is a sum over spin 
and the dispersion relation is given by 



sg: 



n(p3)(P3 + 2|eS|n)i/2 + s/iE 



(27) 



The first term in the square brackets may also be written 
as ^3(1 + 2|ei3|n/p§)^/^ without the sign function. The 
constant an,s ensures that the lowest Landau level only 
contains one spin component, 



f n > 0, 
Ss+ n = 0, eB > 0, 
6s- n^O, eB <0. 



(28) 



We also note again that the phase space associated with 
Landau levels is quantized in a box with periodic bound- 
ary conditions. We omit this to avoid bothersome nota- 
tion like [eL^LyB /2TT \ / LxLy in the phase space factor. 

Let us introduce a constant gauge field A3. One might 
think that a constant gauge field could be gauged away, 
but this is not possible by a gauge transformation satisfy- 
ing the periodic boundary condition. The current density 
is the derivative of the thermodynamic potential with re- 
spect to A3 at the point A3 — 0, 



J3 



dA. 



(29) 



The thermodynamic potential is still given by Eq. (pS)) . 
but the dispersion relation Eq. (j27p is now modified by re- 
placing p3 by Pa -1-6^3. In order to regularize the ultravio- 
let divergences of thermodynamic potential we introduce 
a momentum cutoff A on the P3 integral. Furthermore we 
introduce a cutoff N on the sum over the Landau levels. 
After we have introduced this regularization we can pull 
the derivative with respect to A3 through the sum and 
integral. Then we can use that 



d 
dA3 



dp3 



(30) 



when acting on an arbitrary function of cjp As a result 
we find the following expression for the current density. 



\eB\ ^ ^ dp3 d 



s=± n=0 



27r dp3 



-|-T^log(H-e-'3('^''-±^))j, (31) 



where ujp,s is now given by Eq. (|27|) since we used that 
A3 has to put to zero after taking the derivative. After 
summing over spins the contribution to the integrand of 
the Landau Levels with n > is an odd function of P3. 
Hence only the lowest Landau level which contains one 
spin component contributes to the current. As a result 
for eB > we find 



J3 



\eB\ 



dp3 d 
27r dp3 



OJ, 



,,++T^log(l+e-'^("'''+±^)) 

(32) 



where 



{P3 ± M5) 



(33) 



For eB < one has to replace by _ in Eq. 
Since the integrand is a total derivative, it is easily inte- 
grated. The medium part (logarithmic term) drops be- 
cause it goes to zero with p3 — > ±00. Only a surface term 
remains, which equals 



■^3 ^ K,±(P3 = A) - UJp^±{p3 = -A)] 

\eB\ e^ar 
e^[(A±Ai5)-(ATM5)] = ^B, 



47r2 



27r2 



(34) 



where we have used that ± corresponds to the sign of 
eB. The fact that the current is equal to a surface term 
is because it is caused by the electromagnetic anomaly. 
This as was argued in the first derivation. 

By multiplying the current density Eq. (l34l) with the 
volume one finds the total current Eq. (f24|) . The virtue 
in this derivation is that it is manifest that the current 
results from the surface integral at infinitely large mo- 
mentum, to which any infrared effects of mass, temper- 
ature, and n are irrelevant. The next derivation using 
the derivative expansion will give us more understanding 
why this result is independent of mass. 



D. Derivative expansion of effective action 

The last derivation of the current we discuss is by us- 
ing a derivative expansion of the effective action as is 
performed by D'Hoker and Goldstone [s^l (see also [sl]). 
Let us introduce an axial vector field A^, and write the 
covariant derivative as — dfj, — ieAu 



- ieAl-y^. One 
can define right- and left-handed vector fields as follows: 
Ab. = A^ + A^^ and Al = A^ - A^^. By performing 
the integration over the fermions fields one obtains the 
following effective action 



^eff = logDet {ip 



(35) 



Here Det includes the space-time coordinates as well as 
the color and Dirac indices. The current density can 
be obtained by taking the functional derivative of this 
expression with respect to A^. In the presence of an 
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axial vector field the divergence of a vector current is 
anomalous; one has [531 



F-, 



(36) 



One can write down an expansion of the current in terms 
of the fields A^, A^^ and their derivatives. The expres- 
sion should be Lorentz covariant and U(l) gauge invari- 
ant. Furthermore the current should satisfy the anomaly 
constraint Eq. ([Ml) . To first order in the fields and deriva- 
tives one obtains [53], 



r 



(37) 



The current is m-independent. This follows directly from 
the anomalous divergence of the vector current, that has 
no TO-dependent contributions even with inclusion of a 
mass term. However, the divergence of the axial vector 
current is m-dependent. Therefore the axial vector cur- 
rent induced by a magnetic field depends on mass. This 
is indeed found in Ref. [11]. 

We can now use that eA^ = fj,^ in Eq. ([57]) . so that we 
obtain the current density induced by a magnetic field. 



(38) 



Since the last equation was obtained via a derivative ex- 
pansion, the derivation assumes constant magnetic fields. 



E. Discussion of derivations 

We have argued in Sec. II that A^ ~ d^9/2Nf up to a 
coupling constant. Suppose we have a space-dependent 
theta angle 9, for example formed by a domain wall. The 
covariant current in Eq. ((37|) shows that an electric field 
will induce a current perpendicular to the electric field on 
the domain wall. Moreover, it shows that a magnetic field 
will induce charge on the domain wall. The generation 
of charge on domain walls or solitons was first discussed 
by Goldstone and Wilczek 55]. Callan and Harvey [s^ 
have studied this mechanism as well in the context of 
axionic cosmic strings. They however use pseudoscalar 
coupling instead of axial vector coupling, but find a result 
for the current which is equivalent to Eq. ([57)1 . It was 
argued in Refs. [13, HI] that on domain walls formed 
in certain semi-conductors currents could be generated 
perpendicular to the electric field for the same reason. In 
the context of charge separation in heavy- ion collisions, 
the generation of charge on 9 domain walls was discussed 
by Kharzeev and Zhitnits ky [281] . 

Goldstone and Wilczek [55| have derived their current 
using a perturbative one-loop calculation. It is also possi- 
ble to compute our current perturbatively. One obtains 
a triangle one-loop diagram with two vector couplings 
and one axial vector coupling. As is well known, this 
diagram contains the anomaly. If one includes the effect 



of the chiral chemical potential in the fermion propaga- 
tor, the diagram to compute is the photon polarization 
tensor. 

The axial anomaly generates the topological term 
which is a color singlet. So no net color is separated 
by the Chiral Magnetic Effect. Hence it is expected that 
no additional chromo-electric fields are built up along the 
direction of the magnetic field. Therefore a possible glu- 
onic back-reaction can be neglected. This can also be 
inferred from Eq. (|36|) . since it will not be modified by 
the presence of a gluonic background field. As a result, 
the expression for the current Eq. (|38p is correct even in 
the presence of a time-independent gluonic field. 

If the Chiral Magnetic Effect operates in a heavy- 
ion collision, the current is generated in a finite vol- 
ume. Hence charges are separated, so an electric field 
will be built up along the direction of the magnetic field. 
This could cause a back-reaction. We think that in the 
study for the implications in heavy- ion collisions, this 
back-reaction can be neglected, since the electric field is 
small compared to the magnetic field (it only involves 
a few charges, while the magnetic field is created by all 
charges). Furthermore the electric force is small com- 
pared to the gluonic force. 

We have obtained the current for one fermion with 
charge e. In the quark-gluon plasma there are 3 relevant 
quark flavors, up, down and strange with charges qf = 
2/3e, — l/3e and — l/3e which have A^c = 3 colors. The 
total current will be the sum of the contributions of the 
individual ones, which follow from the previous obtained 
expressions by replacing e with qf, summing over flavors 
and multiplying by the number of colors. This results in 



] 9££ 
L 27r 



(39) 



IV. 



CURRENT EXPRESSED IN CHIRAL 
CHARGE 



As we saw in the previous section, the induced current 
is proportional to ^5. The chiral chemical potential /is 
is a parameter which induces an asymmetry between the 
number density of right- and left-handed fermions 715 = 
nR — hl. Since the asymmetry is conserved by varying 
the magnetic field or the temperature, ^5 will depend on 
the magnetic field, temperature, and chemical potential. 
In this section we will compute the conserved quantity 
ns as a function of /is. We then will express /is in terms 
of ns in order to obtain the dependence of the induced 
current on ns. This allows us to make comparisons of 
the magnitude of the Chiral Magnetic Effect in different 
situations. Moreover, it allows us to relate the current to 
sphaleron dynamics, since the change in Nc, is equal to 
—"^Nf times the winding number of the sphaleron. 

In the computation we present here we will neglect the 
effect of the gluons. At very large temperatures, this is 
correct, since the coupling between gluons and quarks is 
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small due to asymptotic freedom. However, at smaller 
temperatures the relation between /is and could be 
modified by gluonic corrections. A perturbative calcula- 
tion and/or a lattice simulation could give insight in the 
relevance of these corrections. We leave the computation 
of gluonic corrections for future investigation. For QCD, 
the results we obtain at zero temperature are therefore 
unreliable. However, since we keep the discussion gen- 
eral, these results could be of use for a system of non- 
interacting fermions. Again, we will take = Nf — 1 
in the computations. At the end we present the high 
temperature QCD result with multiple flavors. 

In contrast to the computation of current, it is diffi- 
cult to perform the full analytical evaluation of the chi- 
ral charge density. This is because transverse non-zero 
modes have contributions unlike the current which origi- 
nates from zero modes only. We will, therefore, consider 
two simple limits analytically; the weak and strong mag- 
netic field cases in order. Outside these limit we will 
resort to a numerical calculation. 



0.2 




I < r I I I 

2 4 6 8 10 

FIG. 2: Chiral Magnetic conductivity as a function of tem- 
perature and chemical potential. The dashed line is the high- 
temperature/chemical potential approximation. 



A. Weak magnetic field limit 

In the weak magnetic field limit (|ei3| < /is) we can ex- 
pand the current in powers of |e-B|//i|. To leading order 
it is enough to compute the total chiral charge density 
"■5 — nji — riL in the absence of a magnetic field. To 
compute the chiral charge density we first construct the 
thermodynamic potential, 



s=± 



(27r)3 



T^log(l 
± 



-/3(wp_s±At) 



with 



is = iP + s^J■5f + rr? 



(40) 



(41) 



where p = \p\. By differentiating the thermodynamic 
potential with respect to /is one finds the chiral charge 
density which reads 



ns 



7r2 X! 



27r2 



dpp' 



2 + SJ3 



1 - ^ n{ujp^s ± m) 
± 



In the massless limit we find 



ns 



(42) 



(43) 



We can now compute the current in a general magnetic 
field which is assumed to be small compared to /fg. Let 
us define the average magnetic field as 



(44) 



assuming = Ly = = L. Let us for a moment 
assume that $ ^ 27r/e such that to good approximation 



we can ignore the effects of current quantization. We will 
discuss this effect in the next subsection. 

For temperatures and chemical potentials smaller /is 
we find /is « (37r^)^/'^ns^^ such that 



(45) 



Let us define the Chiral Magnetic conductivity ctb as 

J 



OB 



{BY 



(46) 



The Chiral Magnetic conductivity now becomes at zero 
temperature. 



(47) 



For temperatures and/or quark chemical potentials 
larger than /is we find /is ~ 3ns /(T^ -I- /i^/vr^). In that 
case the current yields 



J = 



3e2 



27r2 r2+^2/^2 



N5{B). 



(48) 



The calculation shows that the Chiral Magnetic conduc- 
tivity in the high-temperature limit is 



3e^ 



1 



CTB 



27r2 T2+/i2/7r2 



(49) 



We have displayed the Chiral Magnetic conductivity 
as a function of temperature and chemical potential in 
Fig. [2] using Eq. The figure shows that ctb begins 

from Eq. (|47)) and approaches Eq. (|49p as either T or /i 
grows. 

The reason that the current drops as a function of tem- 
perature and chemical potential is that in both cases 
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1/3 




FIG. 3: Current as a function of flux for T = 0. 



/i5 should take a smaller value for a given 715 through 
Eq. (j43p . This is because a medium at finite T and /i has 
fermion distributions with higher momenta that can take 
part in the chiral charge density. Conversely, a smaller /is 
is sufficient to mimic the effect of a given 77.5 at high tem- 
perature/chemical potential, that means that the effect 
of 715 on systems diminishes by temperature and chemical 
potential. The generated current decreases accordingly 
because it is proportional to /i5. 

Let us briefly discuss current quantization here. If the 
size of the magnetic field is large compared to the area in 
which it is confined the effects of current quantization be- 
come important. For example, consider a magnetic field 
which is constant within a tube with radius R and van- 
ishes outside. The effects of current quantization become 
important if the total flux R^B is comparable to the flux 
quantum 27r/e. 

In Fig. [3] we have displayed the current as a function 
of the flux for T = fj, = 0. Clearly one can see the 
quantization of the current. If becomes an integer 
multiple of 2tt another zero mode is available, which 
result in an increase of the current by an amount of 



B. Homogeneous magnetic field 



Now let us investigate the current in a strong magnetic 
field. We now will take a homogeneous field in which 
we can calculate the induced chiral charge as a function 
of fi^. Again we start from the thermodynamic poten- 
tial which in the presence of a homogeneous background 
magnetic field and a nonzero 115 is given by Eq. (|26p . Dif- 
ferentiating the thermodynamic potential with respect to 



/i5 gives the chiral charge, 
\eB\ ^ ^ 



2tt 



s=± n=0 



dp3 diVp^s 
—X 

1 - ^ n(wp,s ± , (50) 



where 



d/i5 t^p.s 



(51) 



In the massless limit the last equation becomes after 
introducing a cutoff A to regularize the integral 



715 



\eB\n5 



tJ 



27r2 

\eB\ 
2tt 



1 + 2 ^ 



2\eB\n 



'1- 



71=1 
00 ^00 



M5 



EEE 

± s=± n=l 



dp3 dw; 



2tt d/i5 



n{ujp^s ± fJ.)- (52) 



For very large magnetic fields (|ei3| > M5/2) only the 
lowest Landau level (only the first term 1 in the brack- 
ets) contributes to the current. Hence fi^ = 2Ti^n^/\eB\ 
and the current becomes simply equal to the total chiral 
charge in the system. 



J = sgn(B)|e|Ar5 if |eB| > (27r'^ 



)V3„2/3_ 



(53) 



This result can be easily understood from Eq. (Ill|l . In 
a very high magnetic field all modes are fully polarized 
so that we have ^us l) = sgn{eB)nfi L- Applying 
this to Eq. (HJ) gives Eq. ^ [2l|. We shaU limit our 
discussions to the T = ^ = case for a while. If \eB\ < 
/ig/2, not only zeroth but also higher order Landau lev- 
els start to contribute. In Fig. [4] we have displayed the 
current calculated numerically as a function of B. The 
current is saturated for eB/lnzf/^ > {2Tr^y/^ w 5.797. 
The small magnetic field limit result (I45p can be written 



J = 



(Stt^) 



2U/3 



27r2 



-eN. 



eB 
2/3 



(54) 



This limit is displayed as well in Fig d] by a dashed line. 
The approximation is good as long as the current is not 
saturated. The slope is {S-k^)^/^ /{2t:^) « 0.1567. 

The higher-order Landau levels are creating oscilla- 
tions in the conductivity through 715 as can be seen from 
Fig. O where we have displayed the conductivity versus 
ru,^^/\eB\. Related oscillations exist in the conductiv- 
ity induced by an electric field in the presence of a per- 
pendicular magnetic field. In that case they are called 
Shubnikov-de Haas oscillations. The period of the os- 
cillations is equal to ^^5/2 as a function of \eB\. Since 
/i5 depends on \eB\, the period is not a constant func- 
tion of l/\eB\. In the small magnetic field limit we can 
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FIG. 4: Current at zero temperature in a homogeneous mag- 
netic field as a function of magnetic field strength. The dashed 
line indicates the small field limit approximation. 

0.18 I , , , 1 



FIG. 6: Current as a function of magnetic field for different 
temperatures. Displayed with a solid line are from left to 
right: T/ru/^ = 0, 1, 2 and 3. The dashed lines are the small 
field approximations for T jn'^'^ = 2 and 3. 



OB 

^2^2^ri/3 0.16 




0.14 



FIG. 5: Chiral Magnetic conductivity as a function of the in- 
verse magnetic field strength at zero temperature and chem- 
ical potential. 



approximation to be valid for eB < with eB < /i^. 
Indeed this can be seen in the figure, at finite tempera- 
ture the small magnetic field approximation is even good 
to larger values of the magnetic field than at zero tem- 
perature. The oscillations in the Chiral Magnetic con- 
ductivity will be smeared by temperature. 



C. Implications for heavy- ion collisions 

To obtain the induced current in QCD we can sum the 
previous results over flavors and insert the appropriate 
color factor. We will give only the high temperature re- 
sult, since that result will be relatively insensitive to glu- 
onic corrections and is also the most relevant for studying 
the implications of the topological charge changing tran- 
sitions in the quark-gluon plasma. Generalizing Eq. (|48|) 
we obtain 



use that /X5 = (Stt 



2Nl/3„l/3 



Hence for small magnetic 



fields the period of the oscillations in the conductivity 
as a function of nl^'^/\eB\ becomes 2/(37r2)2/3 « 0.2090. 
The mean value of the Chiral Magnetic conductivity can 
be found from the small magnetic field approximation 
which yields ^^/(e^TVs) = {Sn^Y^^ /{2tt^) w 0.1567. 

In order to study the effect of temperature on the cur- 
rent in a homogeneous magnetic field, we have solved 
Eq. (j52p numerically. We have displayed the current in 
Fig. [n]for different temperatures. Clearly, at higher tem- 
peratures it requires larger magnetic fields to saturate 
the current. This is because at high temperature more 
higher momentum modes are occupied, which are more 
difficult to polarize. The dashed lines in Fig.[S]denote the 
small magnetic field approximations from Eq. (j48p . On 
dimensional grounds one expects the small magnetic field 



J 



3e2 Nr, 



1 



27r2 iV/ T2 ^Vtt^ 



(55) 



Here g/e is the electric charge carried by quarks of flavor 
/. Note that Nf in the above is from replacing fi^ by 
7V5. This equation can be used to make predictions for 
the charge asymmetry in heavy-ion collisions like is done 
in Ref. [2^ . In the quark-gluon plasma we could maybe 
expect 77.5 to be several units per fm^ deduced from typ- 

1 /3 

ical QCD sphaleron sizes. In that case T/n^ 1 — 10. 
Since eB ~ 10* MeV^ at the earliest times just after the 
coflision eBjn^^ - 1 - 10; it follows from Fig. ® 
that the current is never expected to be saturated at all. 
Hence the linear approximation in B, Eq. (j55p . can be 
applied to the study of the Chiral Magnetic Effect in 
heavy-ion collisions. 
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In a heavy-ion collision the magnetic field is point- 
ing along the direction of angular momentum, which is 
perpendicular to the reaction plane. We can define like 
in Ref. [2^, to be the difference between the total 
amount of positive/negative charge above and below the 
reaction plane in units of |e|. If /i is small enough we can 
assume that the probability to produce a quark is the 
same as anti-quark. Then each time a sphaleron tran- 
sition with winding number Qw — ^Nc,/2Nf is taking 
place we find that 



/ 



The ± and =F signs in the equations above should be 
read as follows. If the winding number Qw is negative 
(positive), A"*" increases (decreases), while A~ decreases 
(increases). The functions i±{x\_) defined in [2§| are phe- 
nomenological screening functions to describe the effect 
of the quark-gluon plasma through which the separated 
particles have to travel. 

The observables proposed in Ref. [sH and analyzed 
in [H, are sensitive to the correlators (A* A*) 
and (A+A~). These correlators can be obtained from 
Eqs. ([56|) and (|57p by assuming the one-dimensional ran- 
dom walk picture, folding it with the sphaleron rate and 
integrating over time and volume. This analysis has been 
performed in Ref. [1^. 

In Ref. [29I the current was estimated to be propor- 
tional to the degree of polarization of the quarks with 
momenta smaller than the inverse size l/p of the typical 
sphaleron. In Ref. [Soj Eqs. and ([57)1 are similar, 
except that the following replacement has to be made. 



1 



27r2 T2 + /zVtt^ 



(58) 



Since p ~ \/{asT)^ where as is the strong coupling con- 
stant, the newly obtained results are slightly different. 
The difference stems from the fact that in the calculation 
in Ref. poj the typical momenta were determined by the 
inverse size of the sphaleron, a^T, while in the calcula- 
tion in this paper, equilibrium was assumed so that the 
typical momenta are of order T . 



V. EFFECTS OF MASS AND CHIRAL 
CONDENSATE 

In the presence of mass right- and left-handed quarks 
are coupled. A chiral condensate does essentially the 
same. Because the axial charge density operator does 
not commute with the Hamiltonian in the presence of a 
mass term, chirality is not conserved anymore. Hence the 
massive case becomes a dynamical problem and cannot 



be studied using the equilibrium approach we used in this 
article because /is will depend on time when A^s decays. 

The effect of mass on the anomaly was studied in 
Ref. [l^]. It was found that mass term always causes an 
asymmetry between the number of right- and left-handed 
fermions to decay. The time-scale of this decay depends 
the typical momentum (the temperature) of the parti- 
cles, their mass, and the chiral condensate. For T > Tc 
where the momentum scale is much larger than quark 
masses the decay time will be large, so that the equilib- 
rium approach will be reasonably good. However if the 
temperature becomes in the neighborhood of Tc the chi- 
ral condensate becomes important. Any asymmetry will 
be washed out, which will reduce the current. It would be 
very interesting to know how fast the chiral condensate 
washes out the asymmetry; we will leave this problem for 
future study. 



VI. THE CHIRAL BATTERY 

We would like to point out an interesting hypothetical 
application of the Chiral Magnetic Effect - a rechargeable 
battery which stores chirality - the chiral battery. 

Let us imagine a hypothetical material with charged 
fermion quasi-particles described by the massless Dirac 
equation. In this Dirac equation the velocity of light is 
to be replaced by the much smaller Fermi velocity vp oi 
the quasi-particles. A recent example of such a material 
is provided by graphene (for a review see e.g. [s^). even 
though we should keep in mind that chirality in graphene 
is not related to the " usual" spin states considered above 
but instead refers to the sub-lattice states. More directly, 
our considerations may apply to zero-gap semiconductors 
with the linear dispersion relation - possibly, tellurides. 

If we have some finite amount of this material, it can 
be used as a battery. The battery can be charged using 
the axial anomaly by placing it in parallel electric and 
magnetic fields. The charging time will be determined 
by the axial anomaly. The battery stores energy, since 
the Fermi-levels of right- and left-handed modes differ. 
In a sense, this material is also to be regarded as "chiral 
capacitor." 

In the absence of electric and magnetic fields, chirality 
is conserved, so the battery does not discharge. Now let 
us connect the battery to a circuit element with resistance 
R. If we apply a magnetic field to the battery in the right 
direction, a current J will be induced due to the Chiral 
Magnetic Effect. Note that the magnetic field alone does 
no work on fermions in the battery. The behavior of this 
current as a function of the applied magnetic field and 
temperature will follow from our analysis in Sec. HI. The 
current will cause a potential difference V — JR over the 
circuit element. As a result, the same potential difference 
will also exist over the battery. Hence an electric field 
will arise parallel to the magnetic field. In this case the 
axial anomaly operates again to decrease the chirality. 
Hence the rate of discharge will be determined by the 
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axial anomaly as well. 

Let us estimate the amount of energy E stored in 
the chiral battery per unit volume. It is equal to the 
Helmholtz free energy, which is the energy that can be 
used to do work. The free energy is the difference be- 
tween the thermodynamic potential with a chiral charge 
density and without and is easily found by integrating 
Eq. (^5)) with respect to 115. 

We then obtain 



1 



127r 



^4 + l^,l(T^ + !l). (59) 



6' 



At zero temperature and chemical potential we can use 
Eq. to express E in terms of the chiral charge density. 
We find 



E : 



(Stt 



2U/3 



,4/3 



(60) 



The typical distance between the lattice sites in a crystal 
is of order 0.1 nm. Suppose we can store 1 unit of chi- 
rality per lattice site, i.e. an excess of 100 right-handed 
fermions over left-handed fermions per nm^. Then the 
energy density will be 



ij = 7.1xl04^^ 



1.1 X 10 



7 vf J 
c cm'^ 



(61) 



Here is the Fermi velocity. In typical materials like 
graphene vp/c ^ 10~^, so the typical storage capacity of 
the chiral battery is of order 10^ J/cm^ ~ 30 Wh/cm'^. 
This is comparable or better than "conventional" bat- 
teries whose energy density is typically 10 — 100 Wh/Kg; 
note besides that the current in our case is spin-polarized 
and so may be used for spintronic applications. 



VII. CONCLUSIONS 

A system with a nonzero chirality responds to a mag- 
netic field by inducing a current along the magnetic field. 
This is the Chiral Magnetic Effect. The behavior of the 
current as a function of chirality, baryon chemical poten- 
tial and temperature has been obtained in equilibrium in 
this article. 

The Chiral Magnetic Effect can be studied using 
heavy-ion collisions. The possible experimental obser- 
vation of the Chiral Magnetic Effect would be direct ev- 
idence for the existence and relevance of gluon config- 
urations with non-trivial topology. Furthermore it will 
signal V- and CT'-violation in QCD on an event-by-event 
basis. A thorough theoretical understanding of the Chi- 
ral Magnetic Effect will help the experimental analysis 
by offering the possibility of more accurate predictions of 
the observables. 



Since the Chiral Magnetic Effect is due to a mixture of 
QCD and electromagnetic effects, it has very character- 
istic behavior. For example it is expected that the cor- 
relators analyzed in experiment are proportional to the 
square of the charge of the colliding nuclei. This very 
specific behavior can be investigated by measuring colli- 
sions of nuclei with the same atomic number but differ- 
ent charge. With better theoretical understanding, more 
predictions could be made. 

The Chiral Magnetic Effect can only operate in the 
deconfined, chirally symmetric phase. Deconfinement is 
necessary, because quarks need to be separated over long 
distances in order for the Chiral Magnetic Effect to work. 
Restoration of chiral symmetry is needed, since a chiral 
condensate always will wash out any difference between 
the number of right- and left-handed quarks. Hence if 
observed the Chiral Magnetic Effect might be used as an 
order parameter for the confinement/deconfinement and 
the chiral symmetry breakdown/restoration transition. 

Because the Chiral Magnetic Effect probes the P- and 
CP-violating interactions in QCD it can help us to get 
a better understanding of the so-called strong CV prob- 
lem. The problem refers to the fact that strong interac- 
tions do not break the V and CV symmetries explicitly 
even though an addition of V- and CP-odd ^-term to 
the QCD Lagrangian is perfectly allowed without spoil- 
ing gauge invariance. The Chiral Magnetic Effect probes 
the configurations which in principle also cause explicit 
V- and CV violation if 9 is non- vanishing. 

The Chiral Magnetic Effect has also a nice analogy in 
the physics of the Early Universe. One mechanism to 
explain the matter-antimatter asymmetry is electroweak 
baryogenesis [l^, [l3| . There electroweak sphalerons in- 
duce via the axial anomaly C- and CP-odd effects. As a 
result baryon plus lepton number is generated. This pro- 
cess is very similar to the Chiral Magnetic Effect. It is 
also quite possible that the Chiral Magnetic Effect itself 
could have an important role in the Early Universe if a 
large magnetic field and/or a non-zero expectation value 
of the axion field were present at that time. 
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